arXiv:hep-th/9808113vl 19 Aug 1998 


THU-98/31 

hep-th/9808113 


COUNTING PLANAR DIAGRAMS 
WITH VARIOUS RESTRIGTIONS 


Gerard’t Hooft 

Institute for Theoretical Physics 
University of Utrecht, Princetonplein 5 
3584 CC Utrecht, the Netherlands 

e-mail: g.thooft@phys.uu.nl 


Abstract 

Explicit expressions are considered for the generating functions concerning the 
number of planar diagrams with given numbers of 3- and 4-point vertices. It is 
observed that planar renormalization theory requires diagrams with restrictions, 
in the sense that one wishes to omit ‘tadpole’ inserions and ‘seagull’ insertions; 
at a later stage also self-energy insertions are to be removed, and hnally also the 
dressed 3-point inserions and the dressed 4-point insertions. Diagrams with such 
restrictions can all be counted exactly. This results in various critical lines in 
the X-g plane, where A and g are effective zero-dimensional coupling constants. 
These lines can be localized exactly. 


1. Introduction 

In Quantum Chromodynamics (QCD), with N color degrees of freedom and gauge coupling 
constant g, the limit N ^ oo and 0, such that N g'^ = g^ is kept hxed, leads to a 
g expansion, for which all feynman diagrams are planar ^. The 1/N expansion, at higher 
orders in A, requires diagrams on planes with non-trivial topology. 

In this paper, we only consider the g expansion at inhnite N. If ^ is sufficiently 
large, this expansion diverges. Since the theory is asymptotically free, g is small in the 
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ultraviolet region, so that, there, small diagrams dominate. An important question is, 
how the infrared region is affected by the fact that the g expansion has a finite region of 
convergence. Naively, one might expect an infrared Landau ghost, but it is more likely that 
a delicate rearrangement mechanism will take place that will strongly affect the physical 
spectrum of states. One might conjecture that this rearrangement mechanism may cause 
quark co nfi nement in the iV —cx) limit.ref2 

To this end, we study the critical effects for large ^’s in a zero-dimensional model. 
The ‘field theory’ is described by the action 

S{M) = TT{-^M‘^ + lgM^ + lXM^), ( 1 . 1 ) 

where M is an x dimensional matrix, in the linit 

A^ —cx) , ( 7 , A —0, N g'^ = g'^ and NX = X fixed. (1-2) 

Henceforth, the tilde (~) will be omitted. The main body of the paper consists of the 
calculation of the generating functions for the number of planar diagrams with various 
types of restrictions on them. The motivation for the restrictions stems from a special 
renormalization program that is particularly suitable for planar QCD. This program is 
deferred to Appendix A, so as not to interrupt the most important theme of the paper. 

We study many different cases, all related one to another by exact equations. It 
is important to have a consistent notation. Unfortunately, existing notations were too 
haphazard, so we had to invent our own. This is explained in Section 2. 

Much of this paper is based on pioneering work by Koplik, Neveu and Nussinov ^ - 
which in turn makes use of earlier work by Tutte ^ - and on the work by Brezin, Itzyk- 
son, Parisi and Zuber ^. The latter apply matrix theory to do the integration with the 
action (1.1). Here, we choose, instead, to work directly from the functional equations, as 
these will be easier to handle in the QCD case, and they are also more transparent in 
diagrammatic approaches. The relations are read off directly from the diagrams. 

A delicate problem then is the choice of boundary conditions for these equations. 
They can be derived by carefully considering the holomorphic structure that the gener¬ 
ating functions are required to have. Once this is understood, a fundamental solution is 
obtained for the generating function describing the numbers of all planar diagrams for all 
multiparticle connected Green functions, with given numbers V 3 of three-point vertices, V 4 
four-point vertices, and E external lines. It is derived in Sect. 3. From this function, all 
other cases can be derived (Sections 4-8). 

Combining the results, we find the regions in the g, X plane where the planar diagram 
summations (over different kinds of planar diagrams), 

Na{Vs,V 4 ,E)g^^X^\ (1.3) 

V3,V4 

converge. Here, Na, a = 0,1, 2,... refers to the number of diagrams with different kinds 
of restrictions, labled by the number a. 
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An accurate picture of the resulting regions is shown in Fig. 7. Conclusions are 
summed up in Sect. 10. 

Many of our results were obtained and/or checked using the computer program Math- 
ematica. 


2. General notation 


We dehne 

Cs = 

^ 3-point vertices, 


V4 = 

^ 4-point vertices, 


E = 

^ external lines. 


( 2 . 1 ) 


For the generating function for the number of Green functions that include the disconnected 
diagrams, we use the symbol G\ 


CXD, OO, OO 

G{ri,g,\)= ^ G(^e,V3,V4) ■ ( 2 - 2 ) 

E = 0, 

V3=0, V4=0 

Most often, we shall concentrate on the generating functions for connected diagrams only. 
They will be denoted by the letter F: 

OO, OO, OO 

F(z, g,X)^ G^e-i,Vs,V4) ■ (2-3) 

£ 1 = 1 , 

V3=0, V4=0 


Here, it is for technical reasons that we start counting at E = 1, which will become clear. 
Note, also, that we use a different symbol z, instead of g. This is also for later convenience. 

The one-particle irreducible (IPI) diagrams will be generated by the function F: 

OO, OO, OO 

T{u, g,X)^ ■ (2-4) 

E=l, 

V3=0, V4=0 

Apart from limiting ourselves to connected or irreducible diagrams, we can also make 
restrictions on the occurrence of insertions within the (reducible or irreducible) regions of 
the diagrams. This we indicate by using a subscript a, taking values between 0 and 5, 
dehned as in the following list. The symbols G , F , F, and the variables g, A, r], z, and u, 
are all replaced, successively, by Ga, Fa, F^, ga, Xa,rjaT Za, and meaning: 

a = 0, {Go, Fq, ,...): no further restrictions. (2-5) 

a = 1, (Gi, Fi, ,...): no tadpoles: - ^ (2-6) 
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2 : 

no tadpoles and no seagulls: 

(2.7) 

3 : 

also no self-energies: -- 

(2.8) 

4 : 

also no dressed 3-vertices: - 

(2.9) 

5 : 

also no dressed 4-vertices: 

(2.10) 


It will turn out to be important to use the subscripts also for the variables g, X,r],z, and 
u. The reader may appreciate that the use of these suffixes is to be preferred rather 
than exploiting the entire Latin and Greek alphabets. The price payed is slightly more 
complicated-looking expressions. Where ambiguities can be avoided, we will suppress some 
of the subscripts. 

In Eq. (1.1), our action is dehned. Notice, that, in contrast with usual conventions, 
we chose the signs in front of the interaction terms to be positive. Therefore, the positive 
values of the coupling constants are the unstable ones. The reason for this unorthodox 
choice is that now the expansions in terms of Qa and \a all carry plus signs. With this 
sign convention, the coefficients directly correspond to the number of diagrams. Further¬ 
more, unlike Ref ®, we chose factors and in front of the coupling terms. This is 
quite conventional. In our case, these are the best combinatorial factors to choose. The 
coefficients of all expansion terms will be integers, and the only multiplicity factors with 
which diagrams are counted, are the ones corresponding to how many ways a diagram can 
be rotated into distinct forms. In Fig. 1, these multiplicity factors for some diagrams are 
indicated. 



Figure 1. 

Multiplicity factors for some diagrams. 


3. The primary equations 

The recursion equations are most easily derived for the generating function F = Fq de¬ 
scribing all connected planar diagrams. Diagrammatically, we have Fig. 2. Since, in 
this entire section, we will be limiting ourselves to the case a = 0 (no further restric¬ 
tions), we temporarily omit the subscripts 0 that indicate this (so F, g, X, ... stand for 
-^0, go, Ao, ...). 

Writing 

F{z, g, X) = F(^Q'j{g, X) + z F(^i'){g, X) -t- z'^ F{ 2 ){ 9 , A) -|- ... (3-1) 
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Figure 2. Equation (3.2) 


(using brackets to distinguish the subscripts of the Taylor expansion from the subscript 
a), the equation of Fig. 2 reads: 


F^z + g[F^ + -—^ 


+ 


A ( + 2F -—^ + 

z 


F - Fjo) - F{i)Z 

z"^ 


+ F(^o)- 


F-F, 


( 0 ) 


(3.2) 


The last term is very important. If one does not take it into account, the resulting equa¬ 
tions, at a later stage, lead to horrendous complications impeding explicit solution* 

The contributions of F(o) and are necessary for removing unwanted contributions 
from Green functions with too few external lines to close the loops in diagrams (c) and 
(e)~(h) in Fig. 2. 

Clearly, one can solve Eq. (3.2) if ^(0) and F(i) are known. The trick to hnd these 
was described by Koplik et al ^, and goes back to Tutte ^. However, we need the result 
for general g and A, and it appears to be difficult to extend their method directly. The 
matrix integration procedure used by Brezin et al ^, for pure g and pure A theory, 
does work, also for our more general case, but we can easily reproduce this general result 
without integrating over matrices. What needs to be done, as Brezin et al did, is hrst to 
concentrate on the general Green functions Gq , instead of the connected ones, Fq . 

The relation between Ga and Fa is simple to derive diagrammatically. Consider a 
(connected or disconnected) diagram G. Draw a circle around it, such that the circle is 
cut into E pieces by the E external lines. Select one of these E pieces as a starting point. 
This breaks the rotational symmetry, and consequently any diagram with E external lines 
will be counted E/S times, where S is the dimension of the rotational symmetry group 
of the diagram (cf Fig. 1). We open the circle at this point, and then hnd the relation 
between F and G diagrammatically. See Fig. 3. In this Figure, the wiggled line represents 


* In principle, Eqs. of the form (3.16)-(3.21) can then still be derived, but they become 
expressions involving polynomials whose degrees are well over 100, containing coefficients that are 
integers with hundreds of decimal places. 
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Figure 3. Recursive equation relating Ga and Fa 

the circle that was opened up. The cross is the arbitrary point at the boundary that was 
selected. 

The resulting equation, which turns out to hold for all a = 0, ..., 5, is 

OO 

G'(h) = ^ + Yl = 1 + rj G{r])F{r] G{r]), g, X) . (3.3) 

Hence, 


E=1 


F{vG{r]),g,X) + 


gG{g) g 


(3.4) 


Noticing, that g G{r]) acts as z variable in the f un ction F, we rewrite this equation by 
making the following identihcations in Eqs. (2.2) and (2.3): 


Ga — ^ + Fa — -, 

Va 

1 _ ir ^ 1 

— ^ a \ 

Va Za 

(valid for all a = 0 ,..., 5). 

Substituting (3.6) into (3.2) gives us the recursive equation for G{g): 

z + {x — -)(—1 + gx + Ax^) — -{ax + /3) = 0, 

Z .2 


(3.5) 

(3.6) 


where 


1^1 
X = - = F + - ; 
g z 

= ^-^0(0) ; P{9: = fi'-^O(O) + ^-^0(1) + ^-f^l3(0) i 

2 ^gG{g ), 


(3.7) 


(3.8) 

(3.9) 
(3.10) 


and all this only holds for a = 0. The symbol x is the same as the one used to describe 
the Eigenvalues of the matrices M , by Brezin et al ^. 

As this equation is only quadratic^in z, it is easy to solve: 

2^; = X — gx"^ — Xx^ ± \Jx‘^{l — gx — Ax^)^ + A{Xx‘^ + {9 — ot)x + /3 — l) . (3.11) 


^ Here, the last term of Eq. (3.2), diagram (h) in Fig. 2, was crucial. 
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The, as yet unkwown, functions a{g, A) and I3{g, A) can now be determined by repeating 
the arguments by Tutte ^ and Koplik et al ^, and Brezin et al ^. We contemplate the 
branch points by looking at the roots of the 6th degree polynomial in x under the square 
root symbol, as X,g,a and {3 become small. At x small compared to 1/g and 1/V^ 
(which means rj large compared to g and V^)), there are only two roots, 

X — xi^2 ~ ±2 + 0{g, \/X). (3-12) 

All other roots must be at much larger values of x, i.e. rj = 0{g, V^). Branch points 
at such small values of rj are inadmissible, however, and therefore they must cancel out 
pairwise. Thus, for the other roots, one must demand 


x^ — X4^ ; x^ — xq . (3.13) 

This means that it must be possible to write Eq. (3.11) as 

2z — X — gx^ — Xx^ ± A(a: — xs){x — X 5 )^y(x — xi)(x — X 2 ) ■ (3-14) 

It can also be seen from Eq. (3.5) that, for large x, the function z, with the appropriate 
sign choice for the square root, goes as l/x, but we do not need this information here. 

In order to bring the ensuing expressions into a slightly simpler form, we rewrite 

xi + X 2 = 2p ] x\X 2 — —A — Aq . (3.15) 


In view of Eq. (3.12), this implies that p and q will be of order g and A. The requirement 
that Eq. (3.11) can be written in the form (3.14), because of Eqs. (3.13), implies two 
constraints, which hx the functions (y.{g, A) and (3{g, A). Explicitly, what one hnds is: 


9 = 


X = 


a. = 


p (3 + 3p^ — 3^^ + 2p^(^ 


/3 


(1 + ( 3')(6 — 3p2 + p'^ + 12q — 3p2 q-\-Qq^) 

2q — 2p^ + 2q^ — p^q 
(1 + ( 3')(6 — 3p2 + + 12 g — 3p2 g + 6 g^) 

p {2 q — 2p^ + 2 — p^q) (3 + p^ + 3 g + p^ g — 3 + p^ — 3 g^) 

(1 + (3')(6 — 3p2 + p^ + 12 g — 3p2 g + Gg^)^ 

12 g — 3p^ + llp"^ — 5p® + p® — 28p‘^q + 30p^ q — 6p^ q + 

+ 32q^ — AQp^q^ + 19 p"^ q"^ + 2Aq^ — 20 p^ q^ + p'^ q"^ — A q^ 


(6 - 3p2 + p4 + 12 Q - 3p2 g + 6 q2)2 

and with Eqs. (3.9): 

p (3 + p^ + 3 g + p^ g — 3 + p^ — 3 g^) 


A = Fo(o) — 


6 


(1 + g) 


B = To(i) — 


3 p2 _|_ p 4 _|_ 12 g _ 3 p2 g _|_ 6 g2 

^36 - 33p2 + 33 p^ - 13p® +p® + 132 g - 105p2gN^ 
+ 85 p^ g — 18 p® g + 168 g^ — 126 p^ g^ + 77 p^ g^ 

_ 6 p6 g2 _|_ 72 g 3 _ 78 p2 g 3 _|_ 31 p 4 g 3 _ p6 g 3 

- 12 g^ - 33p2 q^ + Qp^q^- 12 g^ - 9p2 g^ 


(6 — 3p2 +p4 + 12g — 3p2 g + 6g2)2 


(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 

(3.21) 
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The way to interpret these equations is that p and q should be eliminated, to hud A and 
B as functions of g and A. 

The elimination process induces singularities as branch points in the 4-dimensional 
space of complex X values. These in turn determine the convergence of the combined 
(( 7 , A) expansion. Since all expansion parameters for To(o) and To(i) must be non-negative 
integers, the singularities at real and positive values for g and A are of particular relevance. 
They will further be discussed in Section 9 and Appendix C. 

Two special cases are: 


1 ) g4>^ theory: 


A = 0 ; 


2 2 q( 1 -Sq) 

V = - 

^ 2 + q 


(3.22) 


In this case, direct solution^leads to a single parameter, a = g/p^ replacing p and q: 


2 2ct 2a 

^ ^ (l-a)(l- 2 a) ’ ^ “ l- 2 a ’ 

g^ = ^ = ia(l -a)(l - 2 a) ; 

a(l — 3a) 

13 - g 0 ( 0 ) - 2(1 - 2 a) ■ 


a = 0 ; 


(3.23) 

(3.24) 

(3.25) 


The last two equations determine the function id)(o) (g) , which is here the only function 
needed to solve Eq. (3.2). 


2 ) X(f)^ theory: 


g = 0 ; p = 0 ; a = 0 ; 

(3.26) 

■'“3(1 + 5)^’ 

(3.27) 

a q{3>-q) 

13 - Ar„,„ - , 

(3.28) 


^ Different procedures lead to different auxiliary parameters ^ these can all be related 
one to another by one-to-one mappings. Up to a sign, the parameter a shown here is the one 
used by Brezin et al . 



4. One-particle-irreducible diagrams 


The generating fnnction T for one-particle-irredncible (IPI) diagrams has been treated in 
Ref ^. Here, we briefly review the procedure. The treatments of this Section hold for all 
cases a = 0 ,..., 5, so we suppress the suffix a. 


Let 


r(n) = —u + g+ Xu^ + loop corrections 
= —u + r(o) + f + r('2)'u^ + ..., 


(4.1) 


where r(i) = Tj-i) +1 denotes the self-energy contributions containing at least one vertex. 
The relation between T and F can be seen diagrammatically, see Fig. 4. 



Figure 4. Equation (4.2) for F 

This Figure corresponds to the equation 

F(z) = 2 + r(o| + f, 1 ) F{z) + r, 2 | (^( 2 ))" + ..., 


or: 

-z = r(F( 2 )). 


(4.2) 

(4.3) 


We write 

Pa = -Za = -VaGa; (4.4) 

Ua = Fa = -= - . (4.5) 

Va Za Za 

For a small number of external lines, and if F^o) s-nd F(g) vanish, the relations between 
the coefficients Fp) and Fp) are straightforward: 

—Z = F(1) (F(i) 2 -t- F(2)z‘^ + ^{3)Z^ ■ ■ ■ ) + ^(2) (-^(1)^ + -^(2)^^ ■ ■ ■ ) +^(3) (F(i)Z . . . ) 

(4.6) 

or 

1 ^ Fm ^ ^’(i)-P’(3)-2F,|, 

Td) - -;p^ ; r( 2 , - ^ ; r,3)- 

= ■^[■^4)^(4) +5 F§, - 5F(i)F(2)F(3)] 

^(1) 


F5 

^( 1 ) 


etc. 


(4.7) 
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5. Removing tadpole insertions (a = 1) 


Vacuum expectation values of single field operators often vanish. In QCD, global gauge- 
invariance is sufficient to see that 


{0|.44.W|0) = 0. (5.1) 

In our case, this amplitude is F(0, g, A) = Fo(o) • We now concentrate on the case a = 1, 
where we impose that all contributions of this sort are omitted. Clearly, the last diagram 
(h) of Fig. 2 is now absent. But also diagrams (c) and (g) require modihcations, as they 
might produce the unwanted tadpole contributions. We thus obtain 


-^l(O) — 0 5 

and removing all contributions to F'i(o), we see that Eq. (3.2) is replaced by 


(5.2) 


F[z,g,X} = 2 + + 

+ X + ^ 


F - F(i)Z 


+ 


+ 


F - 


(5.3) 


if 


a = 1. 


Here, F, z, g,X,... all stand for Fi, ^i, Ai,... 

There are two ways to proceed with this equation. First, one can again write xi — 
Fi -|- 1 /^ 1 , to obtain a quadratic equation for zi as a function of xi: 


zi + (xi -)(AiXi 

-21 


gi^i 




(5.4) 


and we can proceed as in Sect. 3. 

It is instructive, however, to hnd the direct relations between Fi, zi, gi,... on the 
one hand, and Fq, zq, go,... on the other: 


-^1 — Qi{Fo — Fo(o)) ; 

, hence 

Qi 

go + 3AoFo(o) 


= 


9i = 


Ql 


Xl Qi{xq Fq(’q)) , 


Ai 


Qt ■ 


(5.5) 

(5.6) 

(5.7) 


The form of the function Qi{g,X) can be read off from the requirement that, after the 
transformation, the Lagrangian again has a kinetic part normalized to one: 

C = -^xl - Igoxl - IXqXq = C + Dxi - ^xj - ^gixl - ^Xixj , (5.8) 
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from which one reads off Eqs. (5.5)-(5.7), and 

QI — 1 — 2goFo(^o) — 3Ao-Fo(o) ■ (5.9) 

The coefficient D is adjnsted so as to obey Eq. (5.2), and the constant C is irrelevant. 

Indeed, one checks that with the snbstitntions (5.5)~(5.7) and (5.9), the Eq. (3.2) 
turns into Eq. (5.3). In practice, this is the safe method: to adjust the transformations in 
such a way that the new recursion equation holds. 

Since the coupling constants gi and Ai are now new functions of p and q, their 
critical values will also be displaced. The fact that the critical values also follow the 
transformation rules (5.5)~(5.7) is not quite self-evident. But it is true, and related to the 
general feature that dFo(o) and dFo(i) vanish as soon as dp and dq are such that dgo and 
dAo vanish. We return to this issue in Sect. 9. 


6. Removing seagull insertions (a = 2) 


The recursion relation for diagrams where tadpole and seagull insertions are removed, 
reads: 


Fiz^g^X) — z + g (F^ + 


F - F(i)Z 


) + 


+ XiF^ + 2F 


F - F^i)Z 


+ 


F - F(^i-)Z - F(r,\z'^ 


X2y 


( 6 . 1 ) 


where now F, z,g, X ... stand for F 2 , Z 2 , g 2 , X 2 , ■ ■ ■ 

As the seagulls only affect the dressed propagators, the relation between F 2 and Fq 
is as in Eqs. (5.5)-(5.7), but with a different choice for Q. We hud that Eq. (6.1) is obeyed 
provided that 


F2 

Z2 


92 


Q2 


2 


Q2iFo — To(o)) ; 
zq 

Q2 ’ 

go + 3AoTo(o) 


so X 2 = Q 2 {xo - Fo^o)) 
Xq 


Ai — 


Q 2 


1 “ 2(7oTo(o) — 3AoTq(-q^ — 2AoTo(i) 


( 6 . 2 ) 

(6.3) 

(6.4) 

(6.5) 
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7. Removing self-energy insertions ( a = 3 ) 


After having removed the tadpoles and the seagnlls, we remove the self-energy insertions. 
The recnrsion relation becomes, in diagrams, Fig. 5. The equation for the case a = 3 
becomes 


F ^ z + + x(^F^ + =^ + 


2F‘^ F - z- z‘^F. 


( 2 ) 


) - CF, (7.1) 


with 


^( 1 ) - 1 


or 


F = z + F, 2 )z'^ + O(z^) 


(7.2) 



Figure 5. Removing self-energies. 


In (7.1), the last term removes everything that might have produced a self-energy 
insertion up front. This includes the seagull contribution, which is therefore removed 
automatically, so that the corresponding term in Eq. (6.1), being —2A FF(i) , could be left 
out. It just readjusts the 2 :-independent number C. 

C is now hxed by the requirement (7.2), 73(1) = 1. Comparing Eq. (7.1) with (5.3) 
hves us the relations 

73 = Q^Fi ; zj, — zxjQz \ F C = Xj ; (7.3) 

91, = SilQ'i '1 A3 = Ai/(53^; (7.4) 

Q 3 ' = iMd) ; {QiQif = i/7^o(i) ■ (7.5) 

It is easy to see that these are merely held renormalizations that replace the dressed 
propagators by bare ones. 
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8. Removing dressed 3-vertices and dressed 4-vertices 


The last two steps are remarkably easy. Recursion relations for the last two cases are not 
so easy to write down, but they are not needed. The amplitudes considered remain the 
same, but only the dehnitions of the coupling constants change. All one does is recombine 
all contributions from subgraphs that express dressed 3-point vertex insertions and use 
these to redehne the coupling constant g. After this, we do the same with the 4-vertex 
contributions to redehne A. See Fig. 6. 




Figure 6. Collecting 3- and 4-point subgraphs. 


This allows us to write: 

F4 = Fa ; Z4 — zs ; (8.1) 

94 = i"3(2) = Q3'i"l(2) = {QlQsfFo^2) = ' ^^.2) 

Ao(i) 

^4 = As , (8.3) 

and when we also remove the dressed 4-vertices: 



— A's ; ^5 — zs ; 

(8.4) 

95 

= 94 = -^o(2)/-Po(i)^^^ ; 

(8.5) 

As 

T. p np 2 ^1(1)^1(3) - 2Fi(2)^ 

= l3(l) = 1^3(3) 2^3(2) = 3 

Fi(i) 



-P()(l)-^0(3) — 2 Fo(2)^ 

7-1 3 

(8.6) 


^o(l) 
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9. Critical couplings 


At high orders, the number of diagrams can be read off from the values of g and A where 
the amplitudes develop their hrst singularities. These occur when the elimination process 
of p and q in Eqs. (3.16)-(3.19) develops branch points. A branch point is a set of values 
for p and q where a small variation (dp, dg) produces a vanishing variation (dg, dA). Or, 


det 


dg dg 
dp dq 
dX ^ 
dp dq 


0 . 


The solution of this equation was found to be 

p = t^; 12-4t^+t* 

9 — -s — 1; 6 — 6t — 3t^ + 2f^ 


(9.1) 


(9,2) 


At these points, we have {ga, Xa) — (fi'a? -^a)? with 

^ _ ±2 (3 - t)t^ (6 - 6t - 3t2 + 2t3)i/2 

“ (12-4t3+t4)3/2 

(2 + 2t-f)(6-6t-3f + 2t^) 

= - ( 12-413 + ( 4)2 - ■ 


(9.3) 

(9.4) 


The most important region, where g > 0 and A > 0, is from t = 1 — y/S (where 
Aq = 0 , (7o = 1/a/ 12v^) to t = 0 (where A^ = 1/12, = 0 ). 

An important non-trivial Ending is that for the critical line near the origin, where 
Eq. (9.1) holds, one also has, for the functions a and (3 defined in Eq. (3.9), 

da jda djd id{3 dg /dg dX idX 

dp / dq dp / dq dp/ dq dp/ dq’’ 

so that, at a critical point where (dp, dg) / 0, but (dg, dA) = 0, one also has (da, d/3) = 0, 
and since all other coupling parameters ga, Xa are built from a and /3, we may conclude 
that the critical values of these other coupling constants directly follow from the critical 
values of the parameters p and q. In other words, Eqs. (9.2) apply to all critical couplings 
ga, Xa, a = 0,..., 5. Although this is what one would have expected on physical grounds 
(the divergence of the diagrams should not depend on the choice of coupling parameters), 
the mathematical reason for this phenomenon is not totally clear to the author. 

The results for all critical lines are deferred to Appendix C. We depict the critical 
lines in Fig. 7. We see the regions of convergence grow as more constraints are imposed 
on the diagrams. The most surprising feature perhaps is the kinks in the curves at g = 0 
(since diagrams assemble in even or odd functions of g, the picture is symmetric). One 
would have expected smooth, tilted lines if A is plotted against g^, not g, so that the 
lines should run horizontally at g = 0 in the g A plane. But such is not the case. The 
complete courses in the (g^. A) plane of the critical curves at all t values is quite complex. 
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Figure 7. Convergence regions in the g, A plane 

The convergence domains for the cases a = 0 to 5 are shown. The picture is 
numerically accurate. Dots of equal colors and sizes correspond to equal values 
of p, q and t. 


showing cusps, among others at = 0. In contrast, only the physical regions shown in 
Fig. 7 appear to be rather featureless, (the hgure is highly accurate). 

At zero g and at zero A, the critical points are simple algebraic numbers, as indicated 
in Fig. 7; only at Ag = 0, the value of g^ is a root of a polynomial of a high degree. 


10. Conclusion 

All interesting sets of planar diagrams can be counted, in the sense that the generating 
functions for their numbers can be found analytically. This provides us with exact de¬ 
scriptions of the critical lines in the g, A plane. These lines show kinks at (7 = 0, where 
one would have expected analytic ( 7 ^ dependence. Fig. 7 shows the convergence regions 
accurately, when g and A are positive and real. The detailed mathematical expressions 
are given in the appendices B and C. 
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Appendix A. Constrncting perturbative planar QCD amplitudes at all orders 
without divergences 

A renormalization scheme that appears to be not so well-known ® can be set np in the 
following way. It works particnlarly well for planar held theories in 4 space-time dimensions, 
bnt it also applies to several other qnantnm held theories, notably QED. Presnmably, QCD 
with Nc — 3 can also be covered along hese lines, bnt some technical details have not been 
worked ont to my knowledge. This appendix gives a brief description. 

Onr scheme consists of hrst collecting all one-particle irredncible 2-, 3- and 4-point 
diagrams, and formally considering the non-local qnartic ehective action generated by these 
diagrams. 

Next, consider all diagrams nsing the Feynman rnles derived from this action, bnt 
with the limitation that only those diagrams that are absolutely ultraviolet vonvergent are 
included. No snperhcially nltraviolet divergent snbgraph is accepted. To be precise, we 
omit all diagrams with 4 or less external lines, as well as all diagrams containing any 
non-trivial snbdiagram with 4 or less external lines. 

Clearly, one expects that the IPI snbgraphs with 4 or less external lines have already 
been taken care of by onr nse of the qnartic ehective action instead of the bare Lagrangian. 
The important issne to address is, whether the connting of all diagrams was done correctly 
so as to onbtain the reqnired physical amplitndes. Bnt this is not dihicnlt to prove: 

Theorem: the above procednre correctly reprodnces the complete amplitndes for the 
original theory. No diagrams are over-cou n ted or nnder-connted. 

The proof of the theorem is by inspection. Take any diagram of the original theory. 
Consider all its 2-, 3-, and 4-point snbgraphs, and ascertain that these can be identihed 
nnambignonsly as contribntions of onr ehective Lagrangian. In practice, one draws circles 
aronnd the snbgraphs as illnstrated in Fig. 6. Over- or nnder-connting can only happen if 
any ambignity wonld arise with the identihcation. Snch ambignities can only be expected if 
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two subgraphs are connected by at least two lines. But if they both would have 4 external 
lines or less, the entire (sub-)diagram would have 4-1-4 —2 x 2 = 4 or less external lines, and 
this implies that this entire combination should itself be counted as one single contribution 
to the effective Lagrangian (see Fig. 8a). 

If we would have tried to continue the procedure by including 5-point irreducible 
subgraphs into a quintic effective Lagrangian (or more), then counting problems would 
arise: 5-1-5 — 2x2>5 (see Fig. 8b) 


a) 




Figure 8. Ambiguity when including a 5-point vertex 


9 


in the effective Lagrangian 



From our theorem, one deduces that 

if all one-particle irreducible 2-, 3-, and f-point vertices are known, all amplitudes can be 
derived without encountering any divergent (sub-)graph(s). 

In a planar theory, the diagrams considered are of our class 5. 

What remains to be done to complete a perturbative computational scheme, is to 
establish an algorithm to compute the irreducible 2-, 3-, and 4-point vertices (and, if they 
occur, the tadpole diagrams as well). Actually, this is simple. Consider a 4-point IPI 
diagram F(pi,p 2 ,P 3 ,P 4 ) ■ Here, p^ are the external momenta, and pi + P 2 + Ps + Pi — 0■ 
Now consider the difference 

F(pi -F k,p2 - k,p3,p4) - F(pi,P2,P3,P4) = kf,A^^{pi,k,p2 - k,p^,p4) . (A.l) 

The underlining refers to the fact that, in the function A^, this external line follows 
distinct Feynman rules. 

If we follow a path inside the diagram, we can consider the entire expression for 
as being built from expressions containing differences. For instance, in the propagators: 

1 _ ^ {-2p - fc)^ 

{p-\-k)^-\-m'^ p'^-\-m3 {{p-\-kY-\-m'^'){p'^-\-m'^) ' 

or else, in the 3-vertices: 

{p + k)u - p,, = k^S^„ . (A.3) 

We notice, that the expressions for A^ are all (superhcially) ultraviolet convergent! Actu¬ 
ally, one may set up unambiguous Feynman rules for A^{pi, k,p 2 , —kps^p^) and observe 
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that this amplitude exactly behaves as a 5-point diagram, hence it is (superhcially) con¬ 
vergent. 

For the 3-point and the 2-point diagrams, one can do exactly the same thing by 
differentiating more than once. In practice, what one finds is, that there is a set of rules 
containing fundamental irreducible 2-, 3- and 4-point vertices, and in addition rules to 
determine their differences at different values of their momenta. The complete procedure 
thus leads to the following situation. 

We start by postulating the so-called ‘primary vertex functions’. These are, not 
only the irreducible 2-point functions § r[ 2 ](p,—p), the irreducible 3-point functions 
r[ 3 ](pi,p 2 , —Pi — P 2 ) and the irreducible 4-point functions r[ 4 ](pi,... ,p 4 ), but also, in 
addition, the difference functions Aj^j(p, fc, —p — k) and Aj^j(pi, fc,p 2 — /c, —pi — P2), and 
finally the functions 17[2], obtained by differentiating A|^j once more: 

^[2\(Pl+(l:P2-q, -Pi-P2)-Aj"2j(pi,p2, -P1-P2) = quUl^^^{pi-q,g,P2, -P1-P2), (-4.4) 

where one of the other external lines, pi or p 2 is underlined. The double underlining is 
here to denote that the two entries are to be treated distinctly (because of the factor , 
the functions 17^2]^ are not symmetric under interchange of k and q). 

These primary vertex functions are derived by first considering the differences for T[ 4 ], 
Apj and at two different sets of external momenta. These expressions are handled as 
if they were irreducible 5-point diagrams. These are expanded in terms of planar diagrams 
where all irreducible subraphs of 4 or less external lines are bundled to form the primary 
vertex functions. At one of the the edges of such a diagram, we then encounter one of the 
functions A^ or . 

This way, we arrive at difference equations for the primary vertices, with at the r.h.s. 
again the primary vertices. The primary vertices r[ 3 ], A|^j and r 2 ] are then obtained by 
integrating A^j, and A^j with respect to the external momenta. This completes the 
procedure to obtain all amplitudes by iteration. 

Technical implementation of our scheme requires that in all diagrams, an unambiguous 
path can be defined from one external line to another. In QED, one may use the paths 
defined by the electron lines. In a planar theory, one may define the paths to run along 
the edges of a diagram. 

Three remarks are of arder: 

1. The procedure is effectively a renormalization group procedure. The functions A^ 
and play the role of beta functions. 

2. The procedure is essentially still perturbative, since the planar diagrams must still be 
summed. Our beta functions are free of ultraviolet divergences, but the summation 
over planar diagrams may well diverge. 

§ Here, the entries in the square brackets [.. .] differ from the ones in the curved brackets (. ..) 
of Sect. (4) by one unit. 
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3. The procedure only works if the integrations do not lead to clashes. This implies that 
it is not to be viewed as a substitute for regularization procedures such as dimensional 
regularization. We still need dimensional regularization if we want to prove that the 
method is unambiguous, which, of course, it is, in the case of planar QCD. 


Appendix B. The power expansions in g and A. 

To illustrate how our expressions count diagrams, we produce the hrst few terms of the 
power expansions in g and A. 

The case a = 0. 

Let us consider the double series expansions of g and A with respect to p and q. Expanding 
Eqs. (3.16) and (3.17), we get 


9o = 


2 - 22 ^ 


7 3 84 11 5 13 ( 

2? +2« -y9 +y9 


p + 


13 


17 




7q^+ 


12 


Q 


29 5 77 . 


p^ + 


(B.l) 


7 25 


29 
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185 


9 tJtJ o 


287 


35 
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p^ + 0{p\ q^) ; 


An — 


1 203 

3 ^- 3 ^ 


44 5 . 

3^+3^^ 


2 gy 


1 o 2 10 3 

■-+g-2g2 + _g3 


28 


3 


1 11 13 o 25 o 85 4 133 . 

-1- q -oH- q^ -gH-g" 

618^ 9^ 9^ 18 ^18^ 


p^+ 


98 


(B.2) 


y 9 ® ) p* + 0 (p®, 9 ^' 


Such expansions will be displayed in a short-hand notation: 


p 

9o{p,q)= P^ 


y \ 9A 24 12 12 24 24 


Ao(pw) = 


1 
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g3 

g^ 

g5 

g6 
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5 

7 

9 

11 

13 \ 
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2 

2 

2 

2 

2 

3 

13 

17 

-7 
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29 

77 

4 

6 

4 

12 

2 

4 

\ 

25 

29 
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24 

12 

12 

24 

24 
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1 

3 

2 

3 

1 

4 

3 

5 

3 

-2 \ 

1 

3 

1 

-2 

10 

3 

-5 

7 

28 
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13 

25 
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(B,3) 


(B.4) 
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Of course, the matrices extend to infinity. We invert the series ^ in x and y: 





Ao 


COO 

X 4 

AE 




^0 

( 2 

18 

180 

1890 

20412 

224532 \ 


t(^o, Ao) = 

^0 

12 

368 

7860 

143424 

2393496 

1 37700640 

(B.5) 


90 

Vl28 

7120 

240768 

6390720 

146382336 3033374832 / 




( 0 

3 

18 

135 

1134 

10206 \ 



9(^0, Ao) = 

9o 

4 

84 

1368 

20196 

283176 

3847176 


(B.6) 


9o 

Uo 

1768 

49656 

1127808 

22584528 

415844280 / 




The first Green functions are then 



9o 1 

fl 

6 

45 

378 

3402 

32076 

\ 


-^o(o)(fi'05 Ao) — 

COO 

4 

92 

1572 

23904 

341928 

4712580 


(B.7) 


yo 

^2 

1424 

40128 

912960 

18297792 

337041648 / 




fl 

2 

9 

54 

378 

2916 \ 



Fo(i)igoi Ao) = 

9o 

3 

48 

630 

7776 

93555 

1111968 


(B.8) 


9o 

^24 

856 

20112 

392040 

6868152 

112295160 / 




The case a = 1. 

From now on, we leave the powers of g and A to be understood. Whether the g's come 
in even or odd powers depends on whether the number of external lines is even or odd. 
Before looking at these numbers, we continue to produce more lists. From Eqs. (5.6) and 
(6.7): 


9i(p.9) 
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-1 
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-1 
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24 


■ 24 

823 

10453 


4183 

25513 


46309 

192 

768 


128 

384 


384 

1 

2 

1 


4 

5 

\ 

3 

3 


3 

3 

\ 

4 

23 

26 


33 

28 


3 

6 

3 


2 


127 

499 

5765 


6883 

28807 / 

36 

36 

144 


72 

144 

/ 


(B,9) 


(B.IO) 


I For the inversion, a few more terms not shown in Eqs. (B.3—4) were needed. Similarly, Equs. 
(B.ll) and (B.12) required a few terms not shown in (B.9,10). The coefficients that we do show 
are the correct first few terms of the infinite series expansions. 
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Inverting this gives: 



( 2 

12 

108 

1080 

11340 

122472 \ 


p(^i,Ai)= j 

6 

143 

2700 

45441 

714528 

10741086 

(B.ll) 

i^35 

1685 

201447 

4 

2425941 
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25670925 497809800 / 
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3 

18 
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q{9iAi) = 1 

' 4 

1 

48 

666 

9072 

120852 

1582416 

(B,12) 


( 16 

592 

14328 

292869 

5416524 

93623769 / 



Now, -Fi(o) = 0. The next Green fnnctions are fonnd by expanding Eq. (5.5), nsing (5.3) 
Or (3.2): 

/I 2 9 54 378 2916 \ 

^i(i)(^i,Ai) = 1 15 189 2268 26730 312741 (S.13) 

\4 132 2925 54432 917973 14535288 J 

/ 1 9 81 756 7290 72171 \ 

^i(2)(^i,Ai) = 4 99 1755 27216 393417 5450733 (S.14) 

\24 1044 28674 635607 12420216 223297074/ 

The case a = 2. Similarly, we fonnd: 

/I 0 1 2 10 42 \ 

F2 (i)(^2 ,A2) = 15 35 228 1540 10439 

\4 60 725 7636 74725 695464 / 

/ 1 3 15 76 420 2409 \ 

^2(2)(^2,A2) = 4 45 435 3818 32025 260799 

\24 540 8370 107877 1245960 13365702/ 


(B,15) 

(B.16) 


The case a = 3: = 1 and 

/l 3 12 55 273 1428 \ 

^3(2)(^3,A3) = 1 15 159 1460 12405 100449 

\3 90 1638 23400 288738 3227490/ 

Concentrating now on the IPI fnnctions: 


^ 3 { 3 ){ 93 , As) = 
r3(4)(fi'3, As) = 
r3(5)(fi'3, As) = 


/O 1 2 6 22 91 \ 

04 44 364 2720 19380 

\1 34 596 7852 88251 896972 / 

/O 0 5 45 315 2040 \ 

05 115 1565 16950 161950 

\1 65 1750 31890 465080 5873405/ 


/O 0 0 2 15 90 \ 

009 168 1938 18240 

\0 6 249 5176 77841 970596/ 


(5.17) 

(5.18) 

(5.19) 

(5.20) 


21 



a = A: r 4 ( 2 ) = 94: cind 


r4(3)(fi'4, A 4 ) — 0 
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4004 (S. 21 ) 

90174 j 

910 \ 

32585 (S.22) 
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90 \ 

7098 (S.23) 
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a = 5; r 5 ( 3 ) = As and 


r5(4)(5'5, As) — 


r5(5)(fi'5, As) — 


/O 0 5 10 25 70 \ 

0 5 15 70 355 1770 

\ 1 5 50 430 3240 22422 ) 

/O 0 0 2 3 6 A 

0 0 9 54 225 882 

\0 6 69 424 2535 14796/ 


(S.24) 


(S.25) 


It is now illustrative to identify the diagrams that are being counted. We see that, 
in our formalism, the counting is efficient: each diagram is counted essentially just with 
the multiplicity factors of Fig. 1. The tadpole diagrams have at most a factor two if they 
differ from their mirror image. See Fig. 9. 


—00 —0 

^ ^ 0 


Figure 9. The six diagrams of To(o) at order QqXq (see Eq. (B.7). 


The two-point Green functions are separated in various classes by our scheme. In 
Fig. 10, we see the 48 diagrams of Fo(i) at order ^r^A, of which 15 belong to ^ 4 ( 4 ), and 
only 5 of those are in F 21 (cf. Eqs. (B. 8 ), (B.13) and (B.15)). 

Similarly, we illustrate the contributions of different diagrams to the irreducible 5- 
point function, to order ^r^A (see Fig. 11). We see from Eqs. (B.19), (B.22) and (B.24) 
that F 3 ( 4 ) has 65 entries, F 4 ( 4 ) has 35, and F 5 ( 4 ) has only 5. The diagrams are shown 
in Fig. 11. Note that the distinction is whether there are non-trivial irreducible 3-point 
subgraphs present, or irreducible 4-point subgraphs. 
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Figure 10. The 48 diagrams contributing to the 2 point function, at order g^X. 
15 of these have no tadpole insertion; and 5 have also no seagull insertion. 



5 


30 


30 


Figure 11. The 64 diagrams of F 3 ( 4 ), at order g^X, of which 35 are in F 4 ( 4 ) and 5 in F 5 ( 4 ). 


Appendix C. The critical lines. 

The critical lines (Fig. 7) are found by substituting Eqs. (9.2) in all our expressions. Here 
we give the outcomes. Starting from Eqs. (3.16)-(3.21), one just plugs in the equations 
(5.7), (5.9), (6.4), (6.5), (7.4), (7.5), (8.3) and (8.6), and uses the values of as they 
follow from Eq. (3.2). We ignore the minus sign in (9.3), which is trivial. The results are 
the following algebraic expressions. 


g^ = 2f{‘X- t){6 -6t- + - 4t^ + ; (C.l) 
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A[) = {2 + 2t - - 6t - 3t^ + 2t^){12 - + t^)-^ . 


iC.2) 


^ _ 2(6 - - 61 - 3+ 213)5/2 

“ 33/2 (144 - 2881 + 192 - 168 - A8t^ + 100 - 8t^ - 17+ A ’ 

(2 + 2t-t2)(6-6t-3t2+ 2^3)5 
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{C.3) 

{C.A) 


92 = 


2t^{6- t 2 ) 2(6 -6t- 3t^ + 2t3)5/2 


A= = 


(336 - 6721 + 384+ 256 - 712 + A8t^ + 308 -56f - 45+ 12’ 

(a5) 

(2 + 21 - ^ 2)(6 - 61 - 3+ 2 13)5 


(336 - 6721 + 384+ 256 - 712 + A8t^ + 308 -56f - 45 + 12 


{C.6) 


^ 2^2 (6 - t2)2('24 _ 72t - I2t2 + -lOt^- lOt^ + 

~ (6-6t-3t2 + 2t3)5 

_ (2 + 2t-t2)(24- 72t-12t2 + 56t3_iot4_iot5 + 3t6)2 

^ ^ (6-6t-3t2 + 2t3)5 • 


{C.7) 

{C.8) 


2t^ {6-t^)^{6-3t^ +t^) 

~ (24- 72t- 12t2 + 56t3 _ iot4 _ 10^5 + 3 ^ 6 ) 3/2 ’ 


(ao) 



(aio) 


9^, = 9l ; (^^. 11 ) 

A^ = (3456 - 31104 1 + 25920+ 55296 - 180576- 44064+ 

+ 247824- 27936- 147672+ 52760+ 38076- 24432- 
- 17621^2 + 4506 ti3 _ 921 ^14 _ 202 + 114 - 18x 

(24 - 72t - 12t2 + 56t^ - lOt^ - lOt® + 3t^)~^ . {C.12) 


The curves of Fig. 7 are obtained by plotting against using t as a parameter. 

In all plots, t runs from 1 — \/3 to 0, except for a = 5, where t runs from —.890145 
to 0. The new lower bound is the closest zero of Eq. (C.12) for A 5 . The physical reason 
why this larger domain is needed is easy to understand: removing the irreducible 4-point 
diagrams that only contain 3-point vertices can only happen by invoking a negative counter 
term in A. 
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